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Abstract. The NA60 experiment has studied low- mass muon pair production in 
proton-nucleus collisions with a system of Be, Cu, In, W, Pb and U targets using a 
400 GeV/c proton beam at the CERN SPS. The mass spectrum is well described by 
the superposition of the two-body and Dalitz decays of the light neutral mesons rj, p, 
lu, rf and </>. A new high-precision measurement of the electromagnetic transition form 
factors of the 77 and to mesons is presented, complemented with a measurement of the 
temperature parameter of the p meson in cold nuclear matter. The pt spectra for the 
w and <p mesons are extracted in the full pr range accessible, up to pt = 2 GeV/c. 
The nuclear dependence of the production cross sections for the rj, oj and <j) mesons 
has been investigated in terms of the power law <t p a oc A", and the a parameter was 
studied as a function of pr- 
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The study of the production of low mass vector and pseudoscalar mesons in p-A collisions 
represents a natural baseline for the heavy-ion observations, providing a reference for 
different observable in cold nuclear matter, as strangeness production as a function 
of A. In addition, the study of the nuclear dependence of particle properties in p-A, 
as the p T spectra and the production cross sections, is an effective tool to understand 
the dynamics of soft hadron interactions. In order to address these items, the NA60 
experiment complemented the In-In data with a high statistics dimuon sample collected 
in p-A, exposing to the beam six target materials: Be, Cu, In, W, Pb and U. In this 
way, a comprehensive and detailed study of the production of the rj, p, u> and (f> mesons 
as a function of A and pr has been performed for the first time. 

A description of the NA60 apparatus can be found for example in [1] . The produced 
dimuons are identified and measured by the muon spectrometer placed after a hadron 
absorber. The latter allows one to select the highly rare dimuon events but induces at 
the same time multiple scattering and energy loss on the muons, thus degrading the mass 
resolution of the measurement made in the spectrometer. To overcome this problem, 
NA60 already measures the muons before the absorber, with a silicon pixel spectrometer. 
The final sample of dimuons retained for the analysis is defined according to the quality 
of the matching between the reconstructed tracks in the two spectrometers: this leads 
to a total sample of ~ 180 000 dimuons. When the identification of the production 
target is requested, however, stricter cuts must be applied, reducing the total sample 
to ~ 80 000 dimuons. The small component of the combinatorial background (jr and K 
decays) is estimated via an event mixing technique and subtracted from the real data, 
with an overall precision of a few percent. More details can be found in [2]. 

The mass spectrum is well described by the superposition of the two-body and 
Dalitz decays of the light neutral mesons 77, p, uj, rj' and </>. Considering the 
target-integrated data sample, the electromagnetic transition form factors of the 77 
and u mesons have been studied through the Dalitz decays r\ — > p + p^j and uj — >■ 
H + H'7i°. These form factors are usually expressed in terms of the pole parametrization 
\F\ 2 = (1 — M 2 /A 2 ) -2 , with the Vector Meson Dominance (VMD) model giving 
theoretical predictions both for A~ 2 and A~ 2 , which have now been tested with the 
NA60 p-A data. While the details of the analysis can be found elsewhere [3], we 
report here the final values, extracted from the fit on the corrected mass spectrum 
shown in Figure 1, currently the most precise measurements available for the A 2 : 
A" 2 = 1.950 ± 0.059 (stat.) ± 0.042 (syst.) and A" 2 = 2.241 ± 0.025 (stat.) 
± 0.028 (syst.), both expressed in (GeV/c 2 ) -2 . They are in perfect agreement with 
the values obtained by the analysis on the NA60 peripheral In-In data [4] and the 
Lepton-G results [9], confirming the agreement with the VMD model for the r) meson 
and the discrepancy in the case of the u meson. An improved description of the u 
form factor, given by an alternative theoretical approach recently appeared in [10], 
still underestimates the part close to the kinematical cut-off. Finally, it was possible 
to isolate and study the line shape of the p meson, which was parametrized with the 
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This parametrization 

is well suited to describe the p line shape in cold nuclear matter, with the Boltzmann- 
like factor (MT) 3 ' 2 e - ^ playing a crucial role in shaping the low mass tail. The T 
parameter is found to be 161 ± 7 (stat.) ±7 (syst.) MeV. This is the first measurement 
of the effective temperature of the p meson in cold nuclear matter: it is in agreement 
with the measurement in peripheral In-In [4], and consistent with the value of 170 MeV 
obtained by statistical model fits on particle ratios in p-p interactions. 
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Figure 1. Fit on the acceptance corrected 
mass spectrum, with the contribution of the 
r] — » /U/X7, uj — > /i/J7r and p — > /i/i processes. 
Details can be found in [3]. 



In order to study the pt spectra, the data sample was divided, target by target, into 
Pt bins of 200 MeV/c. For each target and pt bin, a fit is then performed on the 
raw mass spectrum with the superposition of the expected sources. In this way, a 
raw pt spectrum is extracted for the uj and mesons, which is then corrected for the 
acceptance and efficiency as a function of pt, as given by the MC simulations. The u 
and pt spectra resulting after the correction for the acceptance x efficiency are shown 
in Figure 2 as a function of the transverse mass and p\. The spectra have been 
compared to a thermal-like function diV/(pTdpT) — diV/ (mr dmr) oc exp(— m^/T) . 
One can immediately appreciate how the thermal hypothesis clearly fails in describing 
the whole spectra up to 2 GeV/c, a systematical deviation from the pure exponential 
trend being clearly visible for high m^. On the contrary, the power-law function 
&N/dp\ oc (1 +Pt/pI) i a standard form for a mixture of soft and hard processes, 
is adequate to describe the diV7dp T distribution in the whole range considered in the 
present analysis. No significant trend for the mean value (pr) as function of A can be 
identified. The average pt integrated over the different nuclear targets are: (pt)w = 
0.61 ±0.03 (stat.) ±0.03 (syst.) GeV/c, (p T )</, = 0.70 ±0.04 (stat.) ±0.02 (syst.) GeV/c. 
These results can be compared with the NA27 measurements in p-p collisions at a/s = 
27.5 GeV [5]: (p T )^ A27 = (0.591 ± 0.021) GeV/c and (pt)^ 27 = (0.513 ± 0.030) GeV/c. 
As one can see, the NA27 value for (pt) for the u agrees with the one obtained in the 
present analysis. On the other hand, for the (pt) of the 0, there is a disagreement by 
more than 4 (statistical) standard deviations. Furthermore, the NA60 results clearly 
indicate that (pt)</> > (pt)u, while the NA27 results indicate that (pt)^ < (Pr)u>- 
Concerning the meson, HERA-B measured production in p-A at y's = 41.6 GeV [6]. 
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The pt distributions measured by HERA-B are well described by the power-law function 
cited above, with (pp) values fully compatible with the NA60 ones and no definite trend 
as a function of A. 




Fitting the raw mass spectra target by target allows us to study the nuclear dependence 
of particle production. We first consider the particle ratios integrated over pp, 
normalizing the cross sections to the one of the u meson as shown in Figure 3. 
The Opjou ratio appears to be flat with A and consistent with Opjo^ = 1. The 
ratio averaged over A — indicated by an horizontal line in the figure — is Opjo^ = 
1.00 ±0.04 (stat.) ±0.04 (syst.). This is in agreement with the ratio a p j(y w = 0.98 ±0.08 
measured in p-p collisions at = 27.5 GeV by the NA27 experiment [5]. A rather 
different conclusion can be drawn, on the contrary, for the cr^/a^ ratio, for which the 
p-p measurement by NA27 [5] is also shown: the trend of the data points indicates that 
a very significant strangeness enhancement is already observed going from p-p to p-A 
collisions, with a further increase as a function of A. 
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Figure 3. Production cross section ratios Opju^ and a^/cr^ in full phase space. 
Error bars and shadowed boxes account for statistical and systematic uncertainties, 
respectively. 



The nuclear dependence of the production cross sections has been studied, for each 
meson, normalizing the measured yields to the one of the lightest target, namely the Be 
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one. The resulting trend is satisfactorily described by the power-law ct p a oc A a , leading 
to the following estimation for the a parameter for the cu and <fi mesons, integrated over 
p T : = 0.841 ± 0.014 (stat.) ± 0.030 (syst.), a = 0.906 ± 0.011 (stat.) ± 0.025 (syst.). 
A possible bias affecting these results comes from the fact that the Be target may start 
to deviate from the power law above, because of its low mass number making it more 
similar to an incoherent superposition of single nucleons, rather than to a nuclear system 
having collective properties. Indeed, if one removes the Be target from the analysis, 
the values for the a parameters shift to a w = 0.801 ± 0.037 (stat.) ± 0.041 (syst.) 
and = 0.963 ± 0.026 (stat.) ± 0.017 (syst.), the u moving towards a softer nuclear 
dependence, a harder one being found, on the contrary, for the <fi. 
The available data sample also allowed us to investigate the pi dependence of the a 
parameters. This study has also been performed for the rj meson, for which the pr 
coverage starts from 0.6 GeV/c. The results, reported in Figure 4, clearly indicate 
an increase of the a parameters as a function of p?, which can be related to the so- 
called "Cronin effect" originally observed by Cronin et al. for charged kaons [7]. The 
same observation is reported for the meson by the HERA-B Collaboration [6], whose 
points are in remarkable agreement with the NA60 ones. A similar comparison can 
be established for the t] meson, thanks to the data published by the CERES-TAPS 
Collaboration [8]: here, too, a good agreement is observed between the two sets of data 
points. No comparison was possible for the u, due to the lack of available measurements. 
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Figure 4. pp dependence of the a parameter for the rj, lo and (f> mesons. 
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